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1.  INTRODUCTION 


This  report  contains  the  summary  of  the  scientific  and  engineering  work  performed  as  part  of 
the  development  of  the  High  Energy  Proton  (HEP)  and  Low  Energy  Particle  (LEP)  instruments. 
HEP  and  LEP  are  a  part  of  the  SOBEDS  suite  of  instruments  being  developed  by  Amptek,  Inc. 
under  the  present  contract  (F19628-95-C-0227).  The  purpose  of  the  HEP  instrument  is  to  detect 
incident  high  energy  protons  (15<E<300  MeV)  and  measure  their  energy  spectrum.  The  LEP 
instrument  is  intended  to  measure  the  differential  energy  spectra  of  electrons  (  0.1  <  E  <  1  MeV) 
and  protons  (1  <  E  <  15  MeV). 

The  primary  technical  effort  during  the  first  year  of  the  SOBEDS  contract  has  been  devoted 
to  the  design  of  the  HEP  sensor  head.  The  key  issue  in  this  effort  is  the  choice  of  scintillator 
material  to  be  used  as  the  primary  proton  detector  and  the  size  and  arrangement  of  the  various 
sensor  head  detectors.  Section  2  contains  the  summary  of  the  extensive  experimental  work 
carried  out  to  determine  the  proper  choice  of  HEP  scintillator  material.  The  baseline  HEP  sensor 
head  geometry,  as  well  as  its  expected  performance,  is  described  in  Section  3.  The  baseline 
design  for  the  LEP  sensor  is  described  in  Section  4.  A  summary  of  the  work  and  concluding 
remarks  are  contained  in  Section  5. 
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2.  SCINTILLATOR  RESEARCH 

2.1.  Materials  and  their  Physical  Properties 


The  choice  of  scintillator  material  to  be  used  as  the  primary  HEP  detector  is  critical  to  the 
proper  performance  of  the  instrument.  HEP  is  to  measure  the  differential  energy  spectrum  of  in¬ 
aperture  incident  protons  with  energies  between  15  and  300  MeV  and  the  measurement  is  to  be 
accurate  in  all  parts  of  the  Earth’s  magnetosphere,  including  the  inner  radiation  belt.  In  order  to 
achieve  this  goal  the  HEP  scintillator  must  meet  three  requirements:  1)  high  density  and 
stopping  power,  2)  high  light  output,  and  3)  short  light  pulse  decay  time. 

The  first  requirement  (high  density)  is  necessary  because  dense  scintillator  material  with  a 
high  stopping  power  can  be  used  to  make  a  small  detector.  This,  in  turn,  has  two  benefits:  1)  the 
decreased  volume  of  shielding  around  the  detector  will  result  in  a  smaller  mass  of  the  instrument 
and  2)  the  smaller  sensitive  detector  volume  will  decrease  the  background  counting  rate  while  the 
instrument  is  in  the  inner  radiation  belt.  The  second  requirement  (high  light  output)  is  necessary 
to  achieve  good  resolution  in  the  measurement  of  energy  deposited  in  the  scintillator  by  incident 
protons.  The  third  requirement  (short  decay  time)  is  important  due  to  the  large  event  rate  in  the 
scintillator  in  the  inner  radiation  belt. 


2.2.  Testing  and  Experimental  Work 

We  have  evaluated  a  number  of  scintillator  materials  for  use  in  HEP  (see  Table  1).  Of  the 
materials  listed  in  Table  1,  we  have  obtained  all  but  the  LSO.  The  scintillators  are  in  “flight 
instrument”  size:  2.0  cm  diameter  and  2.0  cm  long.  LSO  is  not  currently  commercially  available 
but  is  expected  to  come  on  the  market  in  the  first  half  of  1997. 

We  have  conducted  extensive  testing  with  y-ray  sources  and  proton  beams  in  order  to 
determine  the  proton  stopping  power,  energy  resolution  and  temperature  dependence  of  the  light 
output  of  the  scintillators.  Proton  stopping  power  measurements  are  not  yet  complete  and  a  final 
series  of  proton  beam  tests  is  planned  for  January  and  March  of  1997.  We  have  also  undertaken 
a  program  of  radiation  hardness  measurements,  to  be  completed  in  early  1997. 

At  this  time,  GSO  and  LuAP  are  the  leading  candidates  for  HEP  scintillators.  Both 
materials  have  high  total  light  output  and  excellent  stability  with  variation  in  temperature.  LuAP 
is  denser  than  GSO  and  has  a  larger  stopping  power  than  GSO,  but  GSO  has  better  energy 
resolution.  However,  LuAP  is  new  material,  still  in  the  development  stage,  and  significant 
improvement  in  its  performance  is  possible. 
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Table  1.  Listing  of  physical  properties  of  studied  scintillators. 


Property 

LuAP 

(LuA103) 

GSO 

(Gd2SiOs) 

PWO 

(PbW04) 

LSO 

(Lu2SiOs) 

BGO 

(Bi4Ge3012) 

Density 

(g/cm3) 

8.4 

6.7 

8.3 

7.4 

7.1 

Light  Output 

(phot./MeV) 

10,90011,100 

7,9001800 

100111 

27,00012,700 

8,0601120 

Availability 

R&D 

Product 

R&D 

R&D 

Product 

Rad.  Hard. 

(rads) 

Not  Reported 

106 

Not  Reported 

105-106 

Resolution 

(@  662keV) 

9% 

11% 

Not  Reported 

13% 

11% 

Background 

(cnts/sec-cm3) 

320 

None 

None 

310 

None 

Hygroscopic 

No 

No 

No 

No 

No 

Light  Output: 

T  dependence 

Slight 

Slight 

Not  Reported 

Not  Reported 

Strong 

Timing: 

T  dependence 

Not  Reported 

Not  Reported 

Not  Reported 

Not  Reported 

Strong 

3.  BASELINE  HEP  INSTRUMENT  DESIGN 


We  have  developed  a  baseline  design  for  the  HEP  instrument.  The  purpose  of  this  design  is 
to  allow  both  theoretical  analysis  of  the  HEP  on-orbit  performance  and  experimental  tests  of  key 
features  of  HEP  operation.  Accordingly,  we  are  in  the  process  of  building  a  HEP  engineering 
model,  based  on  the  baseline  design,  for  testing  with  accelerator  proton  beams.  The  results  of  the 
theoretical  and  experimental  work  with  the  engineering  unit  will  be  used  to  modify  the  baseline 
design  into  the  final,  flight,  HEP  configuration. 


3.1.  Mechanical  Configuration 

The  overall  mechanical  design  concept  of  the  HEP  sensor  consists  of  a  cylindrical  sensor 
head  mounted  either  directly,  or  next  to,  an  electronics  enclosure,  a  4”  x  4”  x  2-3”  A1  box.  In  the 
final  HEP  design  the  location  and  form  factor  of  the  electronics  enclosure  may  be  different, 
although  the  total  volume  should  stay  close  to  40  in3.  The  sensor  head  will  consist  of  the 
detector  assembly  and  Cu  shielding,  in  the  form  of  a  cylinder,  a  front  collimator  and  a  rear  plug. 
The  0.5  cm  thick  Cu  shielding  cylinder,  with  openings  for  the  photomultipliers  (PMT’s)  will 
surround  the  detector  assembly.  The  cylinder  will  shield  the  detectors  from  protons  with 
energies  up  to  60  MeV  protons.  The  front  collimator  and  rear  plug  are  also  made  of  copper  and 
have  a  thickness  of  1.3  cm  (range  of  100  MeV  protons).  The  front  collimator  will  have  an 
aperture  of  0.25  cm2  in  area  with  an  opening  half-angle  of  7°. 

A  scale  drawn  sketch  of  the  HEP  detectors  is  shown  in  Figure  1.  The  opening  half 
angle,  0,  is  7°.  The  sizes  of  the  various  detectors  are  as  follows: 

D1 :  700  pm  thick,  0.25  cm2  area  PIPS  detector 

D2:  700  pm  thick,  0.25  cm2  area  PIPS  detector 

S 1 :  LuAP  or  GSO  scintillator  cylinder,  1 .5  cm  in  diameter,  2  cm  long 

D3 :  700  pm  thick,  1 .50  cm2  area  PIPS  detector 

S2:  LuAP  or  GSO  scintillator  cylinder,  2  cm  in  diameter,  2  cm  long 

D4:  700  pm  thick,  3.00  cm2  area  PIPS  detector 

S3:  Fast  Plastic  scintillator  hollow  cylinder,  ID  =  1 .9  cm,  OD  =  2.5  cm 

Notes  on  Figure  1: 

1)  PMT’s  which  view  the  SI  and  S2  scintillators  are  not  shown. 

2)  The  S3  scintillator  is  not  a  complete  cylinder  because  it  must  allow  the  light  guides  that 
connect  SI  with  the  PMT  to  pass  through  it.  S3  will  be  viewed  by  1  or  2  photodiodes  or 
avalanche  photodiodes. 

3)  Cu  shielding  on  one  side  of  the  figure  has  been  omitted  to  show  the  dimension  lines. 

4)  The  actual  collimator  will  have  an  opening  feathered  to  allow  for  a  7°  opening  half-angle 
and  1 .2  mm  A1  degrader  plug  which  will  stop  protons  with  energies  below  1 5  MeV. 
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The  approximate  mass  of  the  instrument  may  be  computed  from  the  following: 

A)  Sensor  Head 

A.l)  Cu  Cylinder:  ID  =  3  cm,  OD  =  4  cm,  length  =  6.5  cm,  density  =  8.9  g/cm3 
VCu  =  6.5*[tc(2.0)2  -  n(1.5)2]  s  36  cm3 
MCu  =  36  cm3  *  8.9  g/cm3  s  320  g 

A.2)  Cu  collimator  and  rear  plug:  dia.  =  4  cm,  length  =  1 .3  cm,  density  =  8.9  g/cm3 
Vw  =  2*tt(2.0)2*  1.3  =  32.7  cm3 
Mw  =  32.7  cm3  *  8.9  g/cm3  =  290  g 

A.3)  SI  and  S2  scint.:  length  =  2  cm,  dia.  =  0.5  and  1 .0  cm,  density  =  8.3  g/cm3 
VSI2  =  2.0*[7x(1.0)2  +  7t(0.5)2]  =  8  cm3 
MS12  =  8  cm3  *  8.3  g/cm3  s  65  g 

A.4)  PMT’ s  (Thom  9112  with  divider  chain) 

Mpmt  =  2  *  70  g  =  140  g 

A.5)  Misc:  Dl,  D2,  D3,  D4,  S3  +  mounting  hardware 
MMisc  =  100  g 


B!  Electronics  enclosure 

The  enclosure  is  similar  in  size  but  slightly  smaller  than  the  CEASE  instrument.  Total 
mass  should  be  somewhat  smaller  than  CEASE. 

MEle  =  900  g. 


The  total  mass  of  the  HEP  baseline  design,  M^p,  will  be 
MHEp  =  MCll  +  Mw  +  MS12  +  Mpmt  +  MMisc  +MEIe 
=  (320  +  290  +  65  +  140  +  100  +  900)  g 
=  1.8  kg  =  4.01b. 


3.2.  Expected  Performance 

The  HEP  instrument  will  identify  incident  protons  and  determine  their  energy  by  the 
pattern  of  detectors  hit  by  the  incident  particle  and  the  energy  deposited  in  the  detectors, 
primarily  in  the  SI  and  S2  scintillators.  In-aperture  proton  events  will  be  classified  into  four 
general  types: 
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Type  A:  proton  traverses  Dl,  D2  and  stops  in  SI, 

Type  B:  proton  traverses  Dl,  D2,  SI,  D3  and  stops  in  S2, 

Type  C:  proton  traverses  Dl,  D2,  SI,  D3,  S2  and  D4. 

Type  D:  same  as  Type  C  but  incident  proton  energy  >  300  MeV. 

In  addition  to  these  four  types  of  events  there  can  be  many  other  types  caused  by  out-of¬ 
aperture  protons  penetrating  the  HEP  shielding. 

The  various  in-aperture  event  types  are  caused  by  protons  with  different  kinetic  energies. 
The  energy  losses  of  protons  in  the  HEP  solid  state  detectors  and  the  SI  and  S2  scintillators  are 
shown  in  Figure  2  and  Figure  3.  Protons  with  energies  above  the  degrader  threshold  of  15  MeV 
and  below  120  MeV  will  stop  in  SI  and  produce  Type  A  events.  Type  B  events  are  caused  by 
protons  with  energies  above  the  Type  A  maximum  of  120  MeV  and  below  175  MeV.  Finally, 
Type  C  events  are  due  to  protons  with  energies  above  175  MeV.  A  list  of  energy  deposition 
patterns  for  in-aperture  events  is  shown  in  Table  2.  Energy  losses  for  event  Type  C  are  for  a 
maximum  proton  energy  of  300  MeV.  Higher  incident  energies  (Type  D  events)  will  produce 
smaller  energy  depositions  (SI  and  S2  <  40  MeV).  Note  that  both  front  and  rear  entry  protons 
will  be  detected  as  Type  C  and  D  events. 


Figure  2.  Proton  energy  loss  in  HEP  solid  state  detectors  (front  entry). 
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Event  Type 

D1&D2 

Enerj 

D3 

iy  Loss  (A)  in 

D4 

MeV 

SI 

S2 

Total 

Energy 

A 

0.8<A 

0 

0 

A<120 

0 

SI 

B 

0.6<A<0.8 

0.8<A 

0 

56<A<120 

A<120 

S1+S2 

C  (Front) 

0.5<A<0.6 

0.5<A<0.8 

0.5<A 

36<A<56 

40<A<120 

See  Text 

C  (Back) 

0.5<A<0.6 

0.5<A<0.8 

0.5<A 

40<A<120 

36<A<56 

See  Text 

D 

0.5<A 

0.5<A 

0.5<A 

A  <40 

A  <40 

See  Text 
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Note  that  there  are  several  pulse  height  checks  that  can  be  made  to  help  discriminate  against 
bad  events.  For  example,  for  a  good  Type  B  event  if  S 1  >  S2  then  S 1 +S2  <  145  MeV  and  if  S2  > 
SI  then  S1+S2  >  145  MeV.  Similarly,  for  a  good  front  entry  Type  C  event,  the  SI  pulse  height 
must  lie  between  30  and  60  MeV. 

The  most  difficult  HEP  measurement  is  that  of  the  incident  proton  energy  for  a  Type  C 
event  (175  <  E  <  300  MeV).  In  this  case,  the  proton  punches  through  all  the  HEP  detectors  and 
does  not  deposit  its  full  energy  in  the  instrument.  The  incident  energy  must  be  reconstructed 
from  the  pattern  of  energy  deposition  in  SI  and  S2.  The  most  straightforward  way  of  making 
this  measurement  is  to  derive  it  from  the  pulse  height  difference  S2-S1  for  a  valid  Type  C  event. 
The  pulse  height  difference  (S2  -  SI)  as  a  function  of  incident  proton  energy  for  200  to  300  MeV 
front  entry  protons  is  shown  in  Figure  4.  The  pulse  height  difference  for  rear  entry  protons  is  the 
same  in  magnitude  but  opposite  in  sign. 

The  spread  (FWHM)  in  pulse  height  distribution  is  between  5  and  10%  of  the  peak  energy 
for  LuAP.  This  corresponds  to  2-4  MeV  at  the  high  energy  (300  MeV)  end  and  3-6  MeV  at  the 
low  energy  (200  MeV)  end.  Thus,  it  should  be  possible  to  subdivide  the  200  to  300  MeV  energy 
range  into  two  bins  200-240  MeV  and  240-300  MeV.  If  the  SI  and  S2  energy  resolution  is 
close  to,  or  better  than,  5%,  three  energy  bins  may  be  possible. 


200  220  240  260  280  300 

Proton  Energy  (MeV) 


Figure  4.  S1-S2  pulse  height  difference  for  front  entry  protons. 
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The  primary  high  energy  proton  population  is  in  the  inner  radiation  belt  and  the  HEP 
instrument  must  function  properly  in  that  environment.  The  most  important  aspect  of  the 
instrument’s  operation  is  that  the  true  energy  spectrum  of  incident  protons  be  accurately 
reconstructed  from  HEP  data.  This  can  be  accomplished  by  a  combination  of  two  techniques:  1) 
minimizing  the  possibility  of  misidentifying  the  energy  of  an  incident  proton  and  2)  by  collecting 
sufficient  “background”  data  to  allow  the  effects  of  misidentified  particles  to  be  subtracted  out 
from  the  measured  HEP  spectra. 

The  proton  fluxes  used  in  this  analysis  of  HEP  performance  were  obtained  from  the 
AP8  model.  The  maximum  values  of  inner  belt  integral  fluxes  as  a  function  of  threshold  energy 
(regardless  of  altitude)  are  listed  in 

Table  3.  If  the  incident  flux  is  isotropic  then  the  omni-directional  flux,  F(E),  is  simply 
related  to  the  unidirectional  flux,  f(E)  by 


f(E) 


F(E) 

4k 


(1) 


Equation  (1)  gives  a  good  approximation  of  f(E)  in  the  equatorial  regions.  At  mid-latitudes  the 
factor  in  the  denominator  becomes  3rr  and  near  the  mirror  points  the  factor  is  approximately  2 n. 


Table  3.  Maxium  proton  fluxes  in  the  inner  belt  (AP8). 


Energy 

(MeV) 

Integral  Flux 
p/cm2 -sec  (x  103) 

5 

1,610 

10 

343 

20 

63.9 

40 

29.9 

50 

25.8 

60 

22.9 

80 

18.3 

100 

14.7 

150 

8.6 

200 

5.0 

250 

3.1 

300 

1.9 

10 


(2) 


The  in-aperture  count  rate,  CR  is  given  by 
CR  (E)  =  GF-f(E) 

where  GF  is  the  geometric  factor  of  the  instrument.  GF  is  computed  from 

GF  =  Ik  J*'"  Aid  )sin(6>  )  d6  =  n  A„  [  1  -  cos2(6>,  )]  (3) 

where  A(0)  is  the,  angle  dependent,  effective  area  of  the  aperture,  A0  is  the  total  area  of  the 

aperture,  and  0m  is  the  maximum  angle  of  travel  which  allows  a  particle  to  reach  the  aperture. 

The  HEP  design  under  analysis  has  A0  =  0.25  cm2  and  0in  =  7°,  hence  GF  =  1 1 .7  x  10‘3  em2-sr. 

Type  A  Event  Counting  Rate: 

CRa  =  GF  *  [f  (E  >  20  MeV)  -  f  (E  >  120  MeV) 

=  (1 1 .7  x  10‘3  cm2-sr)  ( (63.9-1 1 .6)  x  103  p-cm'2-sec"’)  /  (4ti  sr) 

=  49  sec’1 

Type  B  Event  Counting  Rate: 

CRb  =  GF  *  [f  (E  >  120  MeV)  -  f  (E  >  175  MeV) 

=  (1 1.7  x  10'3  cm2-sr)  (  (1 1.6-6. 8)  x  103  p-cm^-sec'1)  /  (47t  sr) 

=  4.5  sec"1 

Type  C  Event  Counting  Rate: 

CRc  =  GF  *  [f  (E  >  175  MeV)  -  f  (E  >  300  MeV) 

=  2  (1 1 .7 x  1  O'3 cm2-sr)  ( (6.8-1 .9) x  1 03 p-cm'2-sec') / (4rt  sr) 

=  9.1  sec'1 

Type  D  Event  (E  >  300  MeV)  Counting  Rate: 

CRd  =  GF  *  [f  (E  >  300  MeV) 

=  2  (1 1.7  x  10'3cm2-sr)  ( 1.9  x  103  p-cm'2-sec'')  /  (4tc  sr) 

=  3.5  sec'1 


Total  Counting  Rate: 

CRt  =  CRa  +  CRB  +  CRc  +  CRd  =  66. 1  sec’1 

The  factor  of  2  for  the  C  and  D  type  events  is  due  to  the  fact  that  both  in-aperture  and  rear  entry 
(180°  from  the  aperture  normal)  high  energy  particles  will  produce  valid  events. 
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3.3.  Background  Counting  Rates 

The  singles  counting  rates  of  the  HEP  detectors  are  the  critical  parameter  in  determining  the 
performance  of  the  instrument.  These  rates,  together  with  the  coincidence  resolving  time, 
determine  the  rate  of  false  event  identification,  which  must  be  kept  well  below  the  true,  in¬ 
aperture,  event  rate.  Since  the  singles  counting  rates  depend  only  on  the  detector  geometry  and 
the  incident  flux,  they  determine  the  required  coincidence  resolving  time. 

3.3.1.  Solid  State  Detectors  (D1,  D2,  D3  and  D4) 

The  geometric  factor  of  a  circular  disk  for  particles  incident  from  the  hemisphere  on  one  side 
of  the  disk  is 

r;r/2 

GF  =  2 7t  Jo  ;40  cos(# )  sin(# )  =  nA0  (4) 

where  A0  is  the  disk  area.  Of  course,  the  disk  has  a  second,  equal  geometric  factor  for  particles 
incident  from  the  other  hemisphere. 

D1  and  D2  Detectors: 

D1  and  D2  are  shielded  on  one  side  by  the  front  collimator  (Eth  =  100  MeV)  and  on  the  back 
side  by  the  SI  scintillator  and  the  side  Cu  shielding.  The  back  side  shielding  is  thick  enough  to 
stop  100  MeV  protons.  The  detectors  are  both  0.25  cm2  in  area. 

CR  =  27tAof(E>100  MeV)  =  2n  (0.25  cm2)  (14.7  x  103  p-cm'2-sec'')  /  (4n  sr) 

=  1,840  sec'1 

The  above  rate  is  for  penetrating  protons  (E  >  100  MeV)  only.  The  D1  singles  rate  is  further 
increased  by  in-aperture  protons  which  penetrate  the  degrader  (15  MeV) 

CR  s  GF  *  f(E>15  MeV)  =  (11.7  x  10'3cm2-sr)  (200  x  103  p-cm'2-sec'')  /  (4iz  sr) 

=  186  sec'1 

Thus  the  D1  and  D2  singles  rates,  CRD12,  are  about  2,000  sec'1. 

D3  Detector: 

D3  is  shielded  on  each  side  by  the  large  scintillators  (SI  and  S2)  so  that  the  minimum  proton 
threshold  energy  is  at  least  100  MeV.  The  detector  area  is  1.5  cm2. 

CRd3  =  27tAof(E>100  MeV)  =  2n  (1.50  cm2)  (14.7  x  103  p-cm2-sec'’)  /  (4%  sr) 

=  11,025  sec'1 

D4  Detector: 

D4  is  shielded  by  the  rear  plug  (Eth  =  100  MeV)  on  the  back  side  and  by  S2  and  the  side  Cu 
shielding  on  the  front  side.  The  front  side  shielding  is  thick  enough  to  stop  100  MeV  protons. 
The  detector  is  3.0  cm2  in  area. 

CRd4  =  2TtAof(E>100  MeV)  =  2n  (3.00  cm2)  (14.7  x  103  p-cm'2-sec'')  /  ( 4n  sr) 

=  22,050  sec'1 
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3.3.2.  Scintillators  (SI,  S2  and  S3) 

The  three  scintillators  have  a  surface  area  in  the  shape  of  a  right  circular  cylinder.  The  flux 
through  such  a  surface  can  be  accurately  estimated  as  follows: 

1)  Consider  a  rectangular  box,  with  a  square  cross  section,  circumscribed  about  the  cylinder. 

2)  The  count  rate  of  particles  through  the  box,  CRB0X,  is 

CRm>x  =  /(£)  2 >4  (5) 

/=I 

where  Aj  is  the  area  of  the  i,h  side.  Note  that  the  single  hemisphere  GF  (jiA0)  for  each  side  is 
used  in  the  calculation.  This  is  to  avoid  double  counting  of  particles  incident  on  the  box. 
Particles  incident  from  the  hemisphere  below  a  given  side  will  strike  one  of  the  other  five  sides 
first.  Thus,  using  the  two  hemisphere  GF  will  double  count  the  particles  traversing  the  box 
because  each  particle  will  be  counted  both  when  it  enters  and  exits  the  box. 

3)  If  the  cylinder  has  a  radius  r  and  length  h,  the  particle  counting  rate  into  the  box  will  be 

CRWX  =  xf(E)[Hr2+rh )]  (6) 

where  the  term  in  the  brackets  is  the  surface  area  of  the  box. 

4)  The  actual  surface  area  of  the  cylinder  is  2n  (  r2  +  rh  ). 

5)  Approximate  the  cylinder  counting  rate  by  CRB0X,  multiplied  by  the  ratio  of  the  surface  areas 


CRcy,  =  C In  /  8)  CRnox  =  0.78^-  f(E)  [8  (  r2  +  r /z )  ]  =  156  F(E)  (r2  +  rh ) 

SI  Scintillator: 

51  has  its  front  face  shielded  by  the  front  collimator  (Eth  =100  MeV),  the  back  face  by  S2 
(Eth  =120  MeV)  and  the  bulk  of  its  body  by  side  Cu  shielding  (Eth  =  60  MeV).  A  small  section, 
underneath  the  light  guide,  will  be  shielded  by  the  light  guide  itself  (Eth  s  40  MeV)  and  by  the  A1 
housing  over  the  PMT’s.  Overall,  assume  that  protons  with  energies  above  60  MeV  will 
penetrate  to  SI  from  all  directions.  The  SI  radius,  r,  is  0.75  cm  and  the  length,  h,  is  2.0  cm. 

CRS1  s  1.56  (22.9  x  103)(0.752  +  0.75  x  2.0)  =  74,000  sec'1 

52  Scintillator: 

S2  has  its  front  face  shielded  by  SI  (Eth  =  120  MeV),  the  back  face  by  the  rear  plug  (Eth  = 
100  MeV)  and  the  bulk  of  its  body  by  side  Cu  shielding  (E,h  =  60  MeV).  A  small  section, 
underneath  the  light  guide,  will  be  shielded  by  the  light  guide  itself  (E,h  =  40  MeV)  and  by  the  A1 
housing  over  the  PMT’s.  Overall,  assume  that  protons  with  energies  above  60  MeV  will 
penetrate  to  S2  from  all  directions.  The  SI  radius,  r,  is  1.0  cm  and  the  length,  h,  is  2.0  cm. 

CRS2  s  1.56  (22.9  x  103)  (1.02  +  1.0  x  2.0)  =  107,000  sec1 
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S3  Scintillator: 


The  S3  scintillator  is  an  annular  detector,  it  does  not  have  front  and  rear  faces.  The  effect  of 
this  geometry  is  to  remove  the  r  term  from  the  brackets  on  the  RHS  of  eq.  (7).  S3  is  shielded  by 
the  side  Cu  shielding  (Eth  =  60  MeV).  The  S3  radius  is  1.25  cm  and  its  length  is  2.0  cm. 

CRS3  s  1.56  (22.9  xlO3)  (1.25x2.0)  =  89,000  sec1 

3.4.  False  Coincidence  Event  Rates 

The  event  rates  for  in-aperture  protons  must  be  compared  to  “false  coincidence”  events. 
In  this  type  of  event,  two  or  more  penetrating  protons  strike  the  HEP  detectors  in  such  a  way  as 
to  be  misidentified  as  a  single  high  energy  proton.  Ideally,  the  false  coincidence  rate  should  be 
negligible  compared  to  the  true  event  rate.  In  cases  where  this  condition  is  not  true,  true  event 
count  rates  and  spectra  may  still  be  recovered  from  the  data  if  1)  the  false  event  rate  is  not 
significantly  greater  than  the  true  event  rate,  and  2)  there  is  a  technique  to  subtract  the  false 
events  from  the  measured  count  rate  and  energy  spectra. 

The  penetrating  protons  are  randomly  distributed  in  time.  The  distribution  of  time  intervals 
between  the  arrivals  of  successive  particles  follows  a  Poisson  distribution  with  the  mean  equal  to 
the  inverse  of  the  incidence  rate.  Under  this  condition,  if  the  coincidence  resolution  time  is  t, 
then  the  false  coincidence  rate  between  m  detectors,  Cf,  is  given  by 

Cj  =  mnxn2...nm  r("'_l)  (8) 

where  n,  is  the  singles  counting  rate  for  the  i1*  detector.  In  the  case  of  two  detectors, 

Cf  -  2  nxn2r  .  (9) 

In  the  subsequent  analysis,  it  will  be  assumed  that  the  coincidence  resolution  time,  x,  is  100 
ns.  The  various  false  coincidence  events  are  summarized  in  Table  4  at  the  end  of  Section  3.4. 


3.4.1 .  Type  A  Events 

In  this  event  type  Dl,  D2  and  SI  fire  and  D3,  S2,  D4  and  S3  do  not.  Expected  true  event  rate 
is  49  sec'1.  Initially,  neglect  the  anti-coincidence  shield  S3  and  consider  only  a  Dl  -  D2-  SI  event. 
Since  Dl  and  D2  are  located  close  together,  assume  that  if  a  particle  strikes  one  of  them  it  will 
also  strike  the  other  one.  In  this  case,  the  random  false  event  rate,  FRA,  is 

FRA1  =  2  (74,000  sec'1)  (2,000  sec'1)  (100  x  10'9  sec)  =  30  sec'1. 

However,  this  is  not  the  only  way  a  false  Type  A  event  may  be  recorded.  Out-of-aperture 
protons,  from  the  full  471  solid  angle,  with  energies  sufficient  to  penetrate  the  shielding,  may 
cause  true  coincidences  as  they  traverse  Dl,  D2  and  SI  on  an  oblique  path.  The  solid  angle  for 
such  events  is  defined  by  Dl  and  the  front  face  of  SI.  The  opening  half-angle  is  about  50°. 
Using  eq.  (3)  to  compute  the  geometric  factor  from  one  hemisphere  yields 

GFd1Si  =  7t  (0.25  cm2)  (1  -cos2(480))  =  0.43  cm2-sr. 
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Since  D1  is  shielded  by  the  front  collimator  (Eth  =  100  MeV)  on  one  side  and  by  SI  and  the  Cu 
shielding  (Eth  =  60  MeV)  on  the  other,  the  expected  false  event  rate  will  be 

FRA2  =  GFD!S1  *  [  f  ( E  >  100  MeV)  +  f  ( E  >  60  MeV)  ]  =  1,290  sec1 . 

It  is  clear  that  the  false  event  rate  is  much  greater  than  the  true  event  rate  and  that  the  anti- 
coincidence  shield  S3  will  be  necessary  if  useful  data  on  protons  with  energies  below  120  MeV 
are  to  be  obtained.  A  veto  efficiency  of  about  95%  will  be  needed  in  order  to  reduce  the  false 
event  rate  to  the  level  of  the  true  event  rate. 

The  true  event  energy  spectrum  can  be  obtained  from  the  measured  energy  spectrum  by 
subtracting  the  false  event  energy  spectrum  (multiplied  by  a  scaling  factor  which  depends  on  the 
D1  and  SI  singles  rates  and  on  the  measured  veto  efficiency).  The  energy  spectrum  of  the  false 
coincidence  events  may  be  obtained  by  recording  the  energy  spectrum  of  Dl-  D2-  SI  -  S3  events. 

False  Veto 

Another  problem  event  is  the  false  veto.  In  this  case,  there  is  valid  Type  A  event  and  a 
random  count  occurs  in  any  one  of  the  other  detectors,  thus  giving  the  signature  of  a  bad  event. 
The  false  veto  rate,  FVA,  is  to  a  good  approximation  the  sum  of  the  coincidence  rates  of  valid 
Type  A  events,  CRA,  and  the  singles  rates  in  the  other  HEP  detectors 

FVa  =  2  CRa  (CRD!2  +  CRD3  +  CRD4  +  CRS2  +  CRS3)  x 

=  2  (50  sec1)  (2,000  +  1 1,025  +  22,050  +  107,000  +  89,000  sec1)  (100  x  10'9  sec) 

=  2.3  sec'1 

This  rate  is  only  about  5%  of  the  true  event  rate. 


3.4.2.  Type  B  Events 

In  this  event  type  Dl,  D2,  SI,  D3  and  S2  fire  and  D4  and  S3  do  not.  Expected  true  event 
rate  is  4.5  sec'1.  The  most  likely  false  event  scenario  is  the  simultaneous  passage  of  two 
particles,  one  through  Dl,  D2  and,  possibly,  SI  and  the  second  one  through  D3  and,  possibly  SI 
and/or  S2.  This  event  rate  is  given  by 

FRb  =  2  CRD12  CRD3  t 

=  2  (2,000  sec'1)  (11,025  sec'1)  (100  x  10'9  sec) 

=  4.4  sec'1 

Note  that  this  calculation  neglects  S3  which  will  veto  out  many  of  these  false  coincidences. 
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SI  Energy  Deposition  (MeV) 

Figure  5.  S1-S2  plane  view  of  proton  energy  deposition  in  SI  and  S2. 

An  additional  discrimination  against  this  type  of  false  event  will  be  provided  by  the  pulse 
height  analysis.  The  allowed  range  of  SI  and  S2  pulse  heights  for  type  B  events  is  very 
restricted  (see  Figure  5)  falling  along  an  almost  straight  line  in  the  S1-S2  plane,  in  particular  SI 
>  60  MeV.  The  pulse  height  distribution  of  the  false  coincidence  event  is  expected  to  be  very 
different.  For  example,  out-of-aperture  particles  which  hit  SI,  D3  and  S2  will  travel  through 
very  oblique  paths  through  the  instrument,  thus  keeping  their  paths  in  SI  and  S2  short.  This  will 
result  in  relatively  small  energy  deposition  in  SI  and  S2,  which  will  place  the  event  in  the  lower 
left  hand  quadrant  of  the  plot  in  Figure  5  and  mark  the  event  as  invalid.  (In  Figure  5,  the  first 
data  point  (0,0)  represents  energy  deposition  by  15  MeV  protons,  each  subsequent  point  is  in  5 
MeV  increment.  Data  points  go  left  to  right  for  Type  A  events,  up  the  Type  B  event  line  and 
down  the  Type  C  and  D  event  lines.  Open  circles  are  for  rear  entry  protons.). 

The  use  of  S3  and  pulse  height  analysis  should  reduce  the  false  coincidence  rate  to  well 
below  the  expected  true  rate  of  4.5  sec'1  and  background  spectrum  subtraction  from  the  measured 
energy  spectrum  should  not  be  necessary.  If  the  background  must  be  subtracted,  the  true  event 
energy  spectrum  can  be  obtained  from  the  measured  spectrum  by  subtracting  the  false  event 
energy  spectrum  (multiplied  by  a  scaling  factor  which  depends  on  the  D1  and  D3  singles  rates 
and  on  the  measured  S3  and  pulse  height  analysis  veto  efficiency).  The  energy  spectrum  of  the 
false  coincidence  events  may  be  obtained  by  a  suitable  combination  of  the  SI  and  S2  energy 
spectra  of  events  with  a  D1  or  D2  singles  trigger  and  events  with  a  D3  singles  trigger. 
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False  Veto 


There  are  two  possible  false  vetoes  for  type  B  events,  random  S3  and/or  D4  counts.  Note 
that  the  effect  of  the  second  veto  is  not  to  get  rid  of  the  event  but  to  erroneously  identify  it  as  a 
Type  C  or  D  event.  The  false  veto  rate  is 

FVb  -  2  CRg  (CRD4  +  CRS3)  t 

=  2  (4.5  sec1)  (22,050  +  89,000  sec1)  (100  x  10'9  sec) 

=  0.1  sec'1 


3.4.3.  Type  C  and  D  Events 

In  this  event  type  Dl,  D2,  SI,  D3,  S2  and  D4  fire  and  S3  does  not.  Expected  true  event  rate 
is  9.1  sec'1  for  C  type  and  3.5  sec'1  for  D  type  events.  Several  different  types  of  false  coincidence 
events  will  be  considered. 

Event  Cl:  Random  triple  coincidence 

Three  independent  particles  strike  Dland  D2,  D3  and  D4  simultaneously.  The  event  rate  is 
FRci  =  3  CRD12  CRD3  CRD4  t 

=  3  (2,000  sec'1)  (1 1,025  sec'1)  (22,050  sec'1)  (100  x  10'9  sec)2 
=  0.015  sec'1 

Event  C2:  Valid  or  false  Type  A  event  in  coincidence  with  rear  entry  proton 

A  valid  or  false  Type  A  event  (see  Section  3.4.1)  occurs  in  random  coincidence  with  a  rear  entry 
proton  which  strikes  D4,  S2  and  D3.  The  solid  angle  subtended  by  D3  and  D4  is  0.43  cm2-sr  and 
the  minimum  energy  needed  to  penetrate  the  rear  shield  and  the  S2  scintillator  is  abut  160  MeV 
so  that 

CRD34  =  GF  *  f  (E  >  160  MeV)  =  (0.43  cm2-sr)  (8.6  x  103  pW-sec'1)  /  (4tc) 

=  295 

The  C2  event  rate  is 

FVC2  ~  2  FRA2  CRD34  t 

=  2  (1,290  sec'1)  (295  sec'1)  (100  x  10'9  sec) 

=  0.08  sec'1 
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Note  that  the  allowed  SI  and  S2  pulse  heights  fall  in  a  small  section  of  the  S1-S2  plane  (see 
Figure  5).  Most  of  the  C2  events  will  fall  outside  the  allowed  range,  thus  permitting  additional 
discrimination  against  these  events. 

Event  C3:  Valid  or  false  Type  B  event  in  coincidence  with  rear  entry  proton 

A  valid  or  false  Type  B  event  (see  Section  3.4.2)  occurs  in  random  coincidence  with  a  rear  entry 
proton  which  strikes  D4  and  S2.  The  event  rate  is 

FVC2  =  2  (CRb  +  FRb)  CRd4  t 

=  2  (4.5  +  4.4  sec'1)  (22,050  sec'1)  (100  x  10"9  sec) 

=  0.04  sec'1 

It  appears  the  false  coincidence  rate  of  type  C  and  D  events  will  be  small  compared  to  the  true 
rate,  thus  allowing  a  determination  of  the  flux  of  high  energy  protons  without  any  background 
subtraction.  The  main  difficulty  with  these  events  will  be  the  reconstruction  of  the  original 
incident  proton  energy. 

3.4.4.  Electron  Induced  Events 

The  A1  absorber  in  the  HEP  aperture  is  0.127  cm  thick,  corresponding  to  the  range  of 
700  keV  electrons.  The  maximum  steady  state,  radiation  belt  electron  integral  flux  above  700 
keV  is  about  6  x  106  e/cm2-sec.  This  corresponds  to  an  in-aperture  count  rate  of 

CRe  =  GF  *  [f  (E  >  0.7  MeV) 

=  (1 1.7  x  10'3  cm2-sr)  (  6.0  x  106  p-cm'2-sec'')  /  (4n  sr) 

=  5,600  sec'1 


Table  4.  Listing  of  HEP  false  coincidence  events. 


Type 

Designation 

Description 

A 

FRa, 

Two  protons  hit  Dl,  D2  and  SI 

A 

FRA2 

Out-of-aperture  single  proton  hits  Dl,  D2  and  SI 

Veto  A 

FVa 

Valid  Type  A  in  coinc.  with  any  other  detector  (false  veto) 

B 

frb 

Two  or  more  protons,  hit  Dl,  D2  and  D3  (also  SI  and  S2) 

Veto  B 

FVb 

Valid  Type  B  in  coinc.  with  any  other  detector  (false  veto) 

C 

FRc, 

Random  triple  coinc.  (D1&D2,  D3  and  D4) 

C 

FRc2 

Valid  or  False  Type  A  event  in  coinc.  with  rear  entry  proton 

c 

FRc3 

Valid  or  False  Type  B  event  in  coinc.  with  rear  entry  proton 
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The  energy  losses  of  energetic  electrons  are  of  the  order  of  270  keV  in  the  solid  state  D1  and 
D2  detectors,  so  that  the  high  threshold  energy  of  250  to  300  keV  should  strongly  discriminate 
against  the  detection  of  electrons.  Thus,  the  inner  belt  electrons  should  pose  no  difficulty  for 
HEP.  Under  extreme  conditions,  following  major  solar  events,  this  rate  could  increase  by  several 
orders  of  magnitude  at  altitudes  above  L  =  1.7  RE.  There  may  be  conditions  under  which  the 
accurate  measurement  of  high  energy  protons  in  the  slot  region  and  the  outer  belt  is  not  possible 
because  the  electron  event  rate  will  be  excessively  high. 

The  Cu  shielding  surrounding  the  instrument  is  sufficient  to  stop  7  MeV  electrons.  Under 
normal  conditions,  the  fluxes  of  these  electrons  are  low  (much  less  than  100  e/cm2-sec)  and  will 
not  cause  any  difficulty.  Under  extreme  conditions,  following  major  solar  events,  the  flux  of 
high  energy  electrons  may  increase  sufficiently  to  cause  measurement  difficulties.  However,  this 
can  only  occur  in  the  slot  and  outer  belt  regions  and  not  in  the  inner  belt. 


3.5.  Electronic  Design 

The  over  all  HEP  signal  processing  scheme  is  as  follows  (see  Figure  6).  The  event  trigger  is 
formed  when  pulses  in  D1  and  D2,  and  possibly  D3,  D4  and  S3  exceed  a  fast  discriminator 
threshold,  set  to  about  250  keV  of  energy  deposition  (4  MeV  for  S3).  The  firing  status  of  these 
detectors  is  converted  to  a  five  bit  pattern,  where  each  detector  is  represented  by  one  bit  (1- 
detector  fired,  0-did  not  fire).  Of  the  32  possible  bit  patterns  only  3  correspond  to  valid  Type  A, 
B,  C  and  D  events.  Other  bit  patterns  represent  “bad”  events  although  some  of  them  may  need  to 
be  measured  in  order  to  determine  the  energy  spectra  and  count  rates  needed  for  background 
subtraction.  The  Coincidence  Register,  CR,  (time  resolution:  100  ns)  determines  if  the  event  bit 
pattern  matches  one  of  the  acceptable  ones  (there  would  be  6  acceptable  patterns,  the  choice  of 
which  patterns  are  acceptable  would  be  under  the  control  of  the  DPU).  If  the  event  bit  pattern 
does  not  match  an  acceptable  pattern,  no  further  processing  is  done.  If  an  acceptable  match  is 
found,  the  Coincidence  Register  signals  the  Event  Trigger  Unit  (ETU)  which  in  turn  1)  sends 
gates  to  the  SI  and  S2  ADC’s,  turning  them  on,  and  2)  sends  a  signal  to  the  Local  Event 
Processing  Logic  (LEPL)  circuitry.  LEPL  performs  a  fast  analysis  on  1)  the  D1-D4  pulse  heights 
(using  the  outputs  of  the  quad  discriminators)  and  2)  the  CR  bit  pattern  and  computes  the  event 
ID  word.  The  event  ID  is  sent  to  the  Local  Lookup  Table  and  Data  Storage  (LLTDS)  where  it 
controls  how  the  digitized  pulse  heights  from  SI  and  S2  are  stored.  Every  5-10  seconds,  the 
DPU  signals  the  LLTDS  and  the  scaler  (which  counts  events  in  each  of  the  seven  detectors)  to 
send  the  accumulated  data,  which  is  further  analyzed  by  the  DPU  and  prepared  for  transmission 
to  the  spacecraft. 

The  current  scope  of  the  Amptek,  Inc.  SOBEDS  effort  does  not  include  the  DPU.  The  DPU 
functions  discussed  above  are  for  a  “generic”  DPU. 
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Figure  6.  Block  diagram  of  HEP  electronics. 
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Figure  7.  Schematic  cross  section  diagram  of  the  LEP  sensor  head. 


4.  BASELINE  LEP  INSTRUMENT  DESIGN 

The  Low  Energy  Particle  (LEP)  instrument  is  intended  to  measure  differential  fluxes  of 
protons  with  energies  between  1  and  15  MeV  and  electrons  with  energies  between  0.1  and  1 
MeV.  LEP  sensor  head  consists  of  four  Si  solid  state  detectors  in  a  telescope  configuration  (see 
Figure  7).  The  Cu  shielding  will  be  approximately  0.5  cm  thick,  enough  to  stop  60  MeV  protons 
and  7  MeV  electrons.  The  A1  foil  at  the  entrance  aperture  will  be  0.009  cm  thick,  preventing 
protons  with  energies  below  800  keV  and  electrons  with  energies  below  30  keV  from  reaching 
the  D1  detector. 

The  upper  part  of  the  LEP  proton  energy  range  will  overlap  with  the  lower  part  of  the  E1EP 
instrument  range.  This  will  allow  on-orbit  cross-calibration  between  these  two  instruments.  It 
may  also  be  possible  to  cross-calibrate  LEP  with  the  DOS  instrument,  although  there  are 
significant  difficulties  with  comparing  data  between  dosimeter  and  telescope  type  instruments. 

Electrons  and  protons  are  differentiated  from  each  other,  and  their  incident  energy  is 
determined  by  the  pattern  of  energy  loss  in  the  Si  detectors.  Proton  energy  losses  in  Si  detectors 
can  be  accurately  calculated  using  analytic  expressions  and  measured  proton  energy  loss  data  (for 
example  Janni’s  tables,  AFWL-TR-65-150).  The  calculated  energy  losses  of  energetic  protons  in 
the  LEP  detectors  are  shown  in  Figure  8.  Protons  with  energies  above  15  MeV  punch  through 
the  D3  and  induce  a  signal  in  the  D4  veto  detector. 
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Electron  energy  losses  in  LEP  detectors  cannot  be  calculated  so  simply  because  of  the  large 
amount  of  scattering  that  electrons  suffer  when  traversing  matter.  We  have  undertaken  a  series 
of  Monte  Carlo  calculations  (Integrated  Tiger  Series  code)  of  the  electron  energy  losses  in  the 
LEP  detectors.  These  calculations  are  being  done  in  collaboration  with  our  consultant,  Dr. 
Stanley  Woolf  of  Arcon  Corp. 

We  are  also  undertaking  a  more  sophisticated  calculation  of  proton  energy  losses  in  the  LEP 
detectors  for  the  same  geometry  as  the  ITS  electron  calculations.  The  results  of  these 
calculations  will  be  used  to  further  refine  the  LEP  design.  A  key  aspect  of  these  calculations  will 
be  to  determine  how  different  sizes  of  the  detector  impact  the  efficiency  and  energy  resolution  of 
the  electron  and  higher  energy  proton  measurements. 

A  series  of  electron  (0.03  to  1.8  MeV)  and  proton  (0.8  to  1.8  MeV)  particle  beam 
measurements  with  a  prototype  LEP  sensor  head  was  carried  out  at  GSFC  in  April  1996.  The 
results  of  this  beam  work  were  used  to  validate  the  LEP  response  calculations  and  test  the 
performance  of  the  sensors.  The  results  of  the  data  analysis  showed  that  the  measured  and 
calculated  LEP  responses  were  in  good  agreement  and  that  the  sensor  performance  was 
satisfactory. 
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5.  SUMMARY  AND  CONCLUSIONS 


The  outcome  of  the  work  performed  during  the  first  year  of  the  contract  has  resulted  in  the 
baseline  designs  for  both  the  SOBEDS  HEP  and  LEP  sensors.  The  HEP  sensor  engineering 
model  is  currently  being  manufactured  for  testing  with  proton  beams.  The  physical  construction 
of  the  LEP  sensor  has  been  deferred  as  a  result  of  PL  request. 

The  HEP  design  has  been  optimized  for  operation  in  the  inner  radiation  belt.  The 
configuration  of  the  detectors  has  been  arranged  to  permit  a  clean  measurement  of  in-aperture 
protons  with  energies  between  15  and  300  MeV  in  the  presence  of  an  intense  flux  of  high  energy 
protons.  This  type  of  measurement  requires  the  use  of  fast,  multiple  detector  coincidence  logic 
and  an  anti-coincidence  scintillator  shield.  The  HEP  engineering  model,  which  embodies  the 
baseline  design  will  be  built  and  tested  during  the  second  year  of  the  SOBEDS  contract.  The 
planned  work  on  the  HEP  engineering  model  will  include  proton  beam  experiments  at  the 
Harvard  Cyclotron  (30  <  E  <  150  MeV)  and  the  Brookhaven  National  Laboratory  Alternating 
Gradient  Synchrotron  (  E>  125  MeV). 

The  LEP  baseline  design  has  been  investigated  using  computer  simulation  methods.  The 
instrument  design  is  sufficiently  straightforward,  so  that  we  have  high  confidence  in  the  accuracy 
of  the  computed  results.  Once  we  are  directed  by  PL  to  build  this  sensor,  we  will  proceed  rapidly 
to  the  manufacture  of  the  LEP  engineering  model  and  perform  the  necessary  electron  and  proton 
beam  calibrations. 
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